Time-resolved optical observation of spin- wave dynamics 
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We have created a nonequilibrium population of antiferromagnetic spin- waves in Cr2 03, and char- 
acterized its dynamics, using frequency- and time-resolved nonlinear optical spectroscopy of the 
exciton-magnon transition. We observe a time-dependent pump-probe line shape, which results 
from excitation induced renormalization of the spin-wave band structure. We present a model that 
reproduces the basic characteristics of the data, in which we postulate the optical nonlinearity to 
be dominated by interactions with long-wavelength spin-waves, and the dynamics to be due to 
spin-wave thermalization. 
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Optical magnetic excitations have been studied exten- 
sively in several magnetic oxides , and have recently 
attracted interest in studies of low-dimensional corre- 
lated electron systems ||f§. ^ n this letter we demon- 
strate how time-resolved, nonlinear optical spectroscopy 
(TR-NLOS) of optical magnetic excitations may be used 
to investigate the interactions and dynamics of short- 
wavelength magnetic modes in strongly correlated sys- 
tems. These short- wavelength excitations are the most 
difficult to treat theoretically, and even in three dimen- 
sions their mutual interaction is not fully understood. We 
present results of femtosecond pump-probe spectroscopy 
of the exciton-magnon [X-M) absorption feature in the 
antiferromagnetic oxide Cr2C>3 ||, with both temporal 
and spectral resolution. This optical absorption feature 
allows us to excite antiferromagnetic spin- waves directly, 
with nonequilibrium occupation distributions weighted 
toward large momenta, high energies, and with suffi- 
cient density to observe interaction effects. We have ob- 
served a novel nonlinear optical effect associated with 
the nonequilibrium occupation dependence of the spin- 
wave (SW) dispersion relation. In semiconductor physics, 
it has been demonstrated that important and nontrivial 
physics of correlated many- particle systems can be clar- 
ified through careful attention to dynamics, using TR- 
NLOS || . Our work indicates that TR-NLOS can be 
used to directly manipulate and study magnetic excita- 
tions in strongly correlated insulators, in addition to the 
charged excitations usually probed with optics. 

In the periodic lattice of a magnetic crystal the con- 
cept of a SW appears naturally when the fcrmionic spin- 
Hamiltonian is expressed in terms of Boson operators 
0. Controlled approximations exist for distributions 
of long-wavelength SWs at low excitation densities, but 
away from these conditions theory must be guided by ex- 
perimental observations. Neutron spectroscopy and lin- 
ear optical spectroscopy have provided some of the best 



evidence for the existence of spin wave renormalization 
(SWR) at elevated temperatures J§|,[|], but as previous 
measurements were largely limited to thermally occupied 
SWs, not much is known about the interactions among 
excitations at short wavelengths. Theory suggests that 
the interactions may undergo qualitative changes as the 
zone boundary is approached ]l0| . In our experiments, 
using laser excitation of the X-M absorption, we are able 
to macroscopically occupy a strongly nonthermal distri- 
bution of SWs, with a distribution heavily weighted to- 
ward the zone boundary. 

The X-M transition may be understood as a magnon 
sideband to an exciton, and is similar in character to 
the well-known two-magnon absorption. In a cubic en- 
vironment, the three electrons per Cr 3+ site possess a 
ground state with 4 A2 symmetry and a lowest excited 2 E 
state at ~ 1.7 eV. In a lattice, these levels couple via su- 
perexchange interactions J and J', of order 50 meV. The 
ground state multiplet develops into antiferromagnetic- 
SWs, and the excited state multiplets into "magnetic ex- 
citon" bands. The total spin projection S z is preserved by 
the two spin excitations, so the spin selection rule Am = 
Amx + Atom = is satisfied. To conserve momentum 
the photon is absorbed by an exciton and a magnon of 
equal and opposite momentum, k x + k M ~ Q. The 
excitation is largely localized to neighboring sites, and 
consequently draws its spectral weight from states over 
the entire Brillouin zone. Neglecting interaction between 
the nearby exciton and magnon for clarity, the X-M ab- 
sorption line shape is given approximately by the joint 
density of states (JDOS) f§: 
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The SW dispersion is stronger than that of the exciton, 
and provides the dominant contribution to the line shape. 
The absorption feature may be qualitatively understood 
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FIG. 1. Pump-probe spectrum at At = 0.3 (heavy 
dash-dotted line), 10 (heavy dotted) and 50 ps (heavy solid), 
with units given on the left, shown together with the absorp- 
tion spectrum (light dashed), with units given on the right, 
and the incident laser spectrum (light solid), shown for refer- 
ence in arbitrary units. 

as the SW DOS, shifted up rigidly in energy by the ex- 
citon energy. This allows us to excite SWs and subse- 
quently monitor their density of states (DOS), using a 
pulsed near-infrared laser. 

We have performed pump-probe spectroscopy of the 
X-M transition, using 100 fs pulses emanating at 76 MHz 
from a Thsapphire laser tuned to Hui w 1.765 eV. The 
beam is split in a 10:1 (pump:probe) ratio and the time 
delay, At = t vro be — tpump is controlled with a delay line. 
Both pump and probe are focused through a microscope 
objective to a 6 /im diameter spot at the sample, which 
is held at 10 K with a cold finger cryostat. We estimate 
the average temperature at the sample to be < 30 K. The 
sample is a 2.2 /im thick film of Cr203 grown epitaxially 
on sapphire, by evaporation of Cr in a reactive oxygen 
environment. The c axis of both Cr203 and sapphire 
is normal to the film surface. At the peak of the X-M 
absorption line, the internal transmissivity of the sample 
is 52%. 

After the sample we measure AT x I inc = I on — I f f , 
the change in the transmitted probe intensity as the 
pump is chopped mechanically. In the small signal 
regime, the normalized change in transmission repro- 
duces the differential absorption of the sample, by the 
relation AT/T « -AaL. We measure AT/T ~ (I on - 
Ioff)/Ioff as a function of wavelength and time delay, 
and substract from this a small background contribution 
which persists at negative time delay because of the finite 
repetition rate. At excitation densities of ~ 10 _3 /Cr, 
we observe the complex spectral feature shown in Fig. |l| 
for three different time delays. The data are well de- 
scribed in terms of two components: a spectrally fea- 
tureless photoinduced absorption (PIA) , which shifts the 
overall AT/T toward negative values and is relatively 
time independent over 100 ps, and a derivative-like line 
shape, which is weakly evident at At = and grows as a 



spectral unit as At increases. 

The magnitude and qualitative line shape shown in 
Fig. [l] may be explained by recognizing that the absorp- 
tion of each photon creates a SW, which in turn renor- 
malizes the overall SW band structure through interac- 
tion. In a cubic ferromagnet with occupation at long 
wavelengths, this renormalization takes the form 

Uk({nv})=ul[l-jj^Y,( nk '} UJ ^ ( 2 ) 

k' 

= - K,E tot ], 

where Stot is the total energy of the excited SWs [Q. 
Similar, more complicated expressions hold for antifer- 
romagnets and for different crystal structures. In our 
experiments, the X-M line shape reflects this renormal- 
ization, as a(u>, {«k}) — Pe-m(^) {"k}) now depends on 
the time dependent distribution of photo-excited SWs, 
and the pump-probe experiment measures the time evo- 
lution of -^-nk(t)- In principle, the exciton disper- 
sion relation should also be renormalized for the same 
reasons, but the exciton bandwidth is a factor of 2 nar- 
rower than that of the SWs, and its dispersion is weakest 
near the zone boundary, so this effect contributes little 
to the overall line shape change. Numerical simulations 
confirm these arguments. Since a is proportional to the 
integrated JDOS given in Eq. (|l]), the total derivative 
includes two distinct contributions, one from the level 
shifts and the other from the change in the inte- 
gration volume associated with the level shifts. Qual- 
itatively, the level shifts produce an overall redshift in 
the energy, while the change in the integration volume 
reduces the spectral weight. 

It is interesting to compare the pump-probe signal at 
long times to the change in the linear absorption induced 
by raising the temperature. The SW energy absorbed 
from the laser is only ~ 2% of the total absorbed energy, 
so from the total incident energy density of 5.6 J/cm 3 
only ~ 100 mJ/cm 3 is absorbed by the SW system. From 
the known SW dispersion relation measured by neutrons 
|l2| , we may calculate the SW contribution to the specific 

heat, C P ~C V = ff ^|T 3 , where D ~ 1.25 x 10" 28 J/cm. 
The energy density of one laser pulse thus corresponds 
to a temperature change of ~ 30 K in the magnetic sys- 
tem. In Fig. ||, we compare the saturated pump-probe 
line shape to that expected from a 30 K temperature 
change, based on the measured temperature dependence 
of the linear absorption. For additional comparison, we 
also show the line shape calculated from equations (|l|) 
and (|^), using a scale factor k derived from the temper- 
ature dependent absorption spectra. The pump-probe 
spectrum has been scaled up by a factor of three, which 
may reasonably be attributed to the simplifications of 

E q : (@). 

The agreement among these curves confirms our as- 
signment of the line shape to SWR, and shows that the 
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FIG. 2. Comparison of the saturated pump-probe line 
shape with simulated and actual changes due to a thermal- 
ized energy density. Solid line: pump probe line shape, x3; 
dotted line: thermal difference spectrum; dashed line: calcu- 
lation using eqs. and (0). 

energy density in the magnetic system at long time de- 
lay is comparable to that absorbed directly by the spin 
system from the pump beam. The magnetic excitons ab- 
sorb the majority of the laser energy, serving as a reser- 
voir. The excess energy is transferred rapidly to quench- 
ing sites, whereupon it is dissipated over 100 ns to sev- 
eral microseconds via nonradiative processes |D| . In the 
steady state, a large number of these defect states will 
be excited, together with the steady-state phonon and 
SW distributions. As a test for indirect energy trans- 
fer to the magnetic system via defects and phonons, we 
have performed two-color pump-probe experiments using 
a Coherent RegA 9000 regenerative amplifier system to 
create high intensity pump pulses at 1.5 eV, well away 
from the X-M absorption feature. We focused a portion 
of this beam onto a 3 mm sapphire crystal to generate a 
white light continuum probe, and used an interference fil- 
ter to select a 15 meV spectral range spanning the X-M 
absorption. For absorbed pulse energy densities ranging 
from 1-100 times those used in the degenerate pump- 
probe experiments, we observed no measureable change 
in the X-M absorption, indicating that the mechanisms 
for energy transfer from defect absorption into the mag- 
netic system occurs on time scales much longer than those 
of interest here. We conclude that the magnetic system 
behaves as a quasi-closed system at least during the first 
nanosecond, and that the dynamics which we observe is 
related to an intrinsic internal thermalization of the op- 
tically induced, nonequilibrium SW population. 

We show the time evolution directly in Fig. |[ where we 
plot the response at different wavelengths, keeping laser 
center wavelength fixed. As in Fig. |l|, we have subtracted 
a small background component which is present at nega- 
tive time delay, which is due to the steady state heating 
of the sample. When the probe frequency is outside the 
A-M-line AT/T exhibits prompt PIA with a decay time 
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FIG. 3. Time-resolved pump-probe signal at different 
wavelengths over energy bandwidth hAcj — 20 meV. Probe 
energy ranges are indicated in the insets accompanying each 
curve, and correspond to the following average energies, listed 
from top to bottom: 1.764 eV, 1.771 eV, 1.754 eV, and 1.759 
eV. 

longer than 500 ps, not shown here. When the probe 
frequency is inside the A-M-line AT/T exhibits both 
prompt PIA and picosecond dynamics. The initial dis- 
tribution of SWs created by the laser are weighted heavily 
toward the zone boundary, because of the factor of 47rfc 2 
in the JDOS integral given in Eq. (|). If all SWs con- 
tributed equally to the renormalization, as suggested by 
Eq. (||), one would expect the SWR and hence the X-M 
line shape to undergo an abrupt change at At = 0, and 
remain unchanged during the internal thermalization of 
the SW system. This clearly is not the case: the initial 
population of k « ir/a SWs contributes little to the X-M 
renormalization line shape. 

The exciton and magnon are initially created on neigh- 
boring sites and should interact moderately with each 
other, but the difference in their relative group velocities 
of ^10 A/ps indicates that they should be well separated 
after a picosecond or less, so we do not believe that the 
observed change is associated with the decay of the X-M 
composite. It is also possible that the process of optical 
absorption in the presence of large k excitations is not de- 
scribed well by the particular approximation used here, 
but must include additional many-particle interactions. 
Such effects, however, would need to suppress the contri- 
bution of SWR to the pump-probe line shape by an order 
of magnitude. The simplest explanation for our result is 
that occupation at large k produces weaker overall SWR 
than those at the zone center. Such strong k dependence 
in SW interactions has long been indicated theoretically, 
and the variation of the interaction at short wavelengths 
may be so strong that the overall interaction effects can- 
cel@. 

Regardless of the underlying reason for the difference 
between zone center and zone boundary SW occupation, 
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we can describe the dynamical response phenomenolog- 
ically by dividing the occupied SW states into two dif- 
ferent populations, those at the zone boundary (6) and 
those at the zone center ( c) , with a boundary in recipro- 
cal space chosen to reproduce the experiments. We as- 
sume that the decay of the initial nonequilibrium energy 
density £/, = X)kek b "k^k into the thermalized energy 
density S c = X^kek n k£k is governed by a single ther- 
malization time, r, and the decay over long times is set 
by an overall energy decay time, T . Clearly, r provides a 
measure of the interactions coupling the zone boundary 
spin- waves to those at lower energy, both directly and via 
phonons. This process may be described by the following 
phcnomcnological rate equations: 

^ = -£»/T-a/r, (3) 

If we further assume that only the energy density due 
to zone center SWs £ c is involved in the renormaliza- 
tion line shape, we obtain the following equation for the 
time-dependent pump-probe signal, valid for each energy 
within the X-M absorption region: 

AT/T = Ol 0(t) + a 2 [l - exp(-t/r)] exp(-t/T), (4) 

where ai and ai are prefactors which depend on the spec- 
tral region of interest. The step function Q(t) is required 
to account for the PIA, which we have taken to be time 
independent, consistent with experiments away from the 
X-M peak. The weak structure at early times, due to 
SWR, is also included in a%, though in principle this may 
be accounted for by an additional term. We have di- 
vided the 10 nm spectral range given by our laser spec- 
trum into ten ranges of equal width, separated by 1 nm, 
and measured the temporal pump-probe response in each 
range. We then found the best global fits to the data of 
Eq. (^), in which a\ and a 2 are allowed to vary with 
probe wavelength, and r = 40 ± 10 ps and T = 2 ± 1 ns 
are constrained to be the same for all ten wavelengths. 
The results for four of them are shown by the theoretical 
curves in Fig. [| This simple model captures very well 
the observed dynamics. 

In summary, we have generated a macroscopic, 
nonequilibrium population of SWs and observed its dy- 
namics, using a pulsed laser spectroscopy. At long times, 
the change in the absorption line shape is well under- 
stood by assuming that the photogenerated SWs induce 
a renormalization of the SW dispersion relation. At short 
times, our results deviate sharply from the predictions 
of this model, indicating that short wavelength SWs do 
not contribute to the renormalization line shape. The 
SWs form a quasi-closed system over 100 ps time scales. 
The evolution of this nonequilibrium population is con- 
sistent with a simple model of SW thermalization, and 



the thermalization time characterizes intrinsic spin-wave 
coupling. The effects reported here provide us with in- 
formation on elementary magnetic excitations that is in- 
accessible through conventional techniques, which typ- 
ically probe only thermally occupied magnetic excita- 
tions. Moreover, the technique may be used quite gener- 
ally in magnetic insulators, and may be applied to a wide 
variety of optical magnetic excitations, including two- 
magnon excitations and the recently assigned phonon- 
bimagnon feature in the undoped cuprates || . By using 
multiple wavelength laser sources, this technique may be 
used in conjunction with excitations across the Mott gap 
to probe directly the interactions between charge and 
spin excitations in magnetic solids. 
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